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Abstract 

We present a derivation of the medium-induced gluon radiation spectrum beyond 
the current limitation of soft gluon emission. Making use of the path integral approach 
to describe the propagation of high-energy particles inside a medium, we study the lim- 
iting case of a hard gluon emission. Analytical and numerical results are presented and 
discussed within the multiple soft scattering approximation. An ansatz interpolating 
between soft and hard gluon emissions is provided. The Landau-Pomeranchuk-Migdal 
QnJ ' effect is observed in the expected kinematic region. 
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1 Introduction 



X 
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In ultrarelativistic heavy-ion collisions, a hot and dense medium is produced. One of 
the most important probes of this medium is the production of hadrons and jets with 



high transverse momentum that undergo energy loss processes. This leads to a collection 



of phenomena, usually referred to as jet quenching (reviews of the main models can be 
found in [IH3])) which have been already observed in heavy-ion collisions at the Relativis- 
tic Heavy Ion Collider (RHIC) (y / J/vjv = 200 GeV), such as the suppression of high-px- 
hadron spectra, and of di-hadron and 7-hadron correlations, see e.g. [TJEJSHS] an( ^ re ^ s - 
therein. 

Now, with the heavy- ion program at the Large Hadron Collider (LHC), the current Pb- 
Pb collision energy is much larger (y / s/v r /v = 2.76 TeV) and, thus, observables characterized 
by semi-hard and hard scales are produced more abundantly. Some striking results have 
already been measured: the observation of a strong asymmetry in dijet transverse energy 
[THS], and a strong nuclear suppression of high-py light [13 [11] and heavy [12] hadrons. 
The magnitude of these features at high energies strongly indicates the importance of a 
good description of high transverse momentum observables like jets. 
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The most commonly used jet quenching models describe the changes in hadron spectra 
and jet properties in terms of radiative energy loss processes |13 [I141[TTh20] . This kind of 
approaches neglects the energy transferred to the medium (recoil) and usually works in 
the limit of soft gluon emissions i.e. calculations of the medium-induced gluon radiation 
spectrum are made in the limiting case where the gluon carries a small fraction x of the en- 
ergy of its parent parton (x — > 0). The implications of such limitation on the computation 
of physical observables are discussed at length in [3J. Nevertheless, the original techniques 
in [13j[l31[T7] can be extended beyond the limit x — > 0. This extension was studied for 
the case of the energy distribution in the BDMPS limi10 in Refs. |14f[IUl[2T?ll21j . Here we 
will present results beyond the BDMPS limit for the case of a hard gluon emission; some 
attempts in a similar direction can be found in [22^23] . 

In this work, we use the path-integral formalism |15Hl~7 ,24-27j to compute a finite-energy 
correction for hard gluon emissions off a quark inside a medium by studying the limiting 
case x — > 1. Afterwards, we generalize the result, in a heuristic manner, to include the 
already known limit of soft gluon emission. The paper is organized as follows. In the next 
section we will provide a small review of the medium-induced gluon radiation spectrum 
in the limit of soft gluon emissions, as well as the formalism used. Our calculation of 
the finite-energy correction will also be made in this section. In section [3] we will present 
the analytic results within the multiple soft scattering approximation |13[I15] for the final 
interpolation expression and the numerical discussion will be made in section [H Finally, 
in section [5] a summary and future prospects will be presented. 

Let us note that information about the size of finite energy corrections to medium- 
induced gluon radiation is badly needed in order to put the existing Monte Carlo models 
for jet quenching, which automatically include energy-momentum conservation - see e.g. 
[28 -30] . on more solid grounds. Together with research on the role of interferences from 
different emitters in the QCD shower [31H35] . and the effects of color exchanges with the 
medium [36, [37], they constitute one of the main avenues in the recent developments of 
the theory of energy loss processes in a QCD medium. 



2 Medium-induced gluon radiation 




q = (1 — x)p 

Figure 1: Radiation process of a gluon off a quark with the corresponding kinematical 
variables. 

Consider the inelastic process shown in figure [TJ A quark of 4-momentuir|^| p, = 0, 

1 We define the BDMPS limit as the one where L — > oo with uj/lo c finite, and lo c — \qL 2 , see [21] 
for details. Here q is the transport coefficient encoding the properties of the medium, L its length and 
lj the energy of the radiated gluon. Formally, this limit is equivalent to integrating the gluon transverse 
momentum up to infinity, neglecting the kinematic constrains. 

2 We will use light-cone coordinates defined as a = (ao, a x , a y , a z ) = (a+, a_ , aj_) with a± — (ao±a z )/\/2 
and a± = (a Xy a y ). Furthermore, Ca = N and Cf ~ (N 2 — 1)/(2N) are the quadratic Casimirs in the 
adjoint and fundamental representations respectively, with N the number of colors. 
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coming from a hard process with amplitude Mh, emits a gluon of 4-momentum k with 
+-component k + = xp+ and transverse momentum kj_. 

Inside a medium, partons undergo multiple scattering. In the high-energy limit, the 
parton essentially conserves its energy and the effect of the medium is only a rotation 
of its color field, thus acquiring an eikonal phase. In this situation, a convenient way of 
describing the parton propagation is using Wilson lines, 



W(xq+, L + ; xj_) = J" exp < ig / dx + A-(x + ,xj_) 

I Jxo+ 



(1) 



which describe the propagation of a particle through a medium with longitudinal bound- 
aries at [cco+, £+] and color field As for the less energetic parton, the restrictions of 
the above formula need to be relaxed to allow some motion in the transverse plane of the 
propagating particle (see e.g. [26|,l27j and refs. therein). By doing this, the Wilson line is 
replaced by the path integral propagator from x + , xqj_ to L + , xj_ , 



G(x 0+ ,x ±;L + ,x_ L \p + ) 



Dr ± (£)exp^-f± 

r±(x +)=^o± [ •'•n,. 

x W(x 0+ ,L + ;r ± (^)), 



dr± 



(2) 



which associates a Brownian motion of the particle in the transverse plane at the same 
time that its color field is rotated. 



2.1 Soft limit: x -»■ 

Considering the soft gluon emission limit (x — > 0), there are two contributions for the total 



amplitude: the case where the gluon does not interact after it is emitted (figure 2(a) ), and 



the case where the radiation vertex is inside the medium (figure 2(b) ). 



q ~ p 
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(a) Radiation vertex outside the medium. 



(b) Radiation vertex inside the medium. 



Figure 2: Radiation diagrams in the limiting case x — > 0. 
This means that the total amplitude is given by 



Mi 



x->0 



tot 



M: 



(3) 



where describes the diagram of figure 2(a) and describes the diagram of 



figure 2(b) Thus, the double differential spectrum of medium-induced radiation in this 
limit reads 



tot 



dl 



rr->0 



\ lvl tot I 
2(2vr) 3 



d jvac 
dk i <i 2 k 



med 



+ k. 



s-»0 



dl 
dk+d 2 \t i 



(4) 
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and makes reference to the 



where (0) denotes the medium average of 0, see e.g. 
spin and color averages. The vacuum contribution is recovered from ^\M^~^°\ 2 y. 

dl vac a s C F 1 



A; 4 



dk+d 2 \t\ 



x^O 



TT 2 k 2 L 



(5) 



The medium contribution comes from the diagram 2(b) with its own complex conjugate 



M% | 2 ) ) , usually called gluon term, and from this same diagram with the complex 



conjugate of diagram 2(a) ({M®~*° (M^ )^ ), denoted interference terrrjl: 

d jmed 



dk+d 2 k 



x-^0 



{2ll) 2 k + \k+ 



dy+dy + dx±e iki ' x± e - 



i/dfn(0«r(x±) 



d 







dy± <9x_l 



X(y+,y± = 0;2/+,x_l|A; 4 



(6) 



+ / dy + dx ± e- ix ^ 2- 







k ± 2 dy ± 

Here, % denotes the two-dimensional path-integral 

rr±(v+)=y± 

%(x + ,x ± ;y + ,y ± \k + ) = / 2>rj_(£) 

Jrj_(x+)=x ± 
[ f V+ 

x exp < / d£ 



X(y+,y±_ = 0;L + ,xJ/c + ) 



(7) 



ik + (dr±\ 1 



This results from the fact that one has to average over all medium configurations since for 
the calculation, only a frozen configuration profile was taken into account. The medium 
dependence of the spectrum comes from the factor n(£)a(r±). The density of scattering 
centers is given by n(£) and their space configuration and strength is contained in the 
dipole cross section, 



a(r ± ) = 2 J 



(2tt)5 



K(q_iJ 



-iq x r 



(8) 



with ao(q±) the scattering potential corresponding to one scattering center, usually taken 
in the form of a static Debye screened potential [38J. 



2.2 Hard limit: x 



We will now to parallel the derivation in the previous section but in the limit x — > 1. This 
means that the kinematics of the process is now constrained by 

|q±l (9) 

We work in the high-energy limit, which means that terms proportional to pi x , qi x or 
k; x are neglected in the numerator of the propagators [26027] . As for the denominators, the 
propagators coming from the initial quark and from the gluon will be simplified to 

inh ith 

(10) 



pi + ie %Pi + Pi_ + ie 



s Full explanations of the variables in this equation can be found in [24H27] . 
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Figure 3: Radiation diagram for gluon emission outside the medium in the one-scattering 
case. p,pi,k,q denote the 4-momenta, x\ the position of the scattering center, and 
A,A±,A' the color indices in the fundamental representation, while a, b denotes the color 
index in the adjoint. 



whereas for the final quark, the term proportional to q^j_ is kept in the denominator: 

id. id 



1 1 



qf + ie 2q i+ qi_ - ctf ± + is 



ill) 



By doing this we are assuming that the most energetic partons (the initial quark and the 
gluon) will acquire only a color phase by crossing the medium. Since the final quark is 
much softer than the other two, the contribution from qj^ in the denominator must be 
taken into account to allow some motion in the transverse plane. 

First, we will compute the contribution where the gluon is emitted outside the medium. 
Starting by one scattering with the medium (diagram represented in figure [3] ) , the 7- 
matrix reads 

(12) 

with A = A b T b . For the simplification of the Dirac structure, the approximations used in 
the soft gluon emission cannot be applied @. As for the integrals, the only non-trivial ones 
are the ones related with the propagators: 



dp- 

e 



- ip .( X(i+1)+ ~ Xi+ ) * = e(a; _ x . +)i ( i3) 



2-7T p_ — (—ie) 

where is the step function. Thus equation (fl2j) can be reduced to 



7f = -gT% Ai |±^± J dx 1+ [ig(A^ AlA (x 1+ ,0 ± )}e(x 1+ )u(q)f^M h (k + q). (14) 

Generalizing to n scatterings, one can check that this structure iterates. Summing over 
all scattering centers an exponential series is found and the T-matrix for a gluon emitted 
outside the medium can be written as follows: 

7 q = -gT% Ax a / *A + W AlA (x 0+ , L+; Q x )u{q)f^M h {k + q). (15) 



4 In the case x — > 0, q ~ p and thus the following simplification was used: u{q)f.*f 1 ~ u{q)f<j[ ~ 2q- e* ~ 
2<jr-|_el, where e — (e+, e_, ej_) is the gluon polarization. Now, for x — > 1, the momenta relation has changed 
and the Dirac equation can no longer be used. Besides, we have to keep j* unevaluated since we do not know 
a priori which component will be dominant as fc+ is larger than in the previous case. This means that the 
Dirac structure can only be reduced to u(q)^*($ + ^)A t A 1 ApMh(p) ~ 2(g + k) + (A-)A 1 AP+u(q)^*'j-Mh(j>) 
(see appendix rXl for the relations between the 7 matrices). 
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Figure 4: Radiation diagram for gluon emission inside the medium for the limiting case 
x — > 1. The meaning of the variables and indices is analogous to that in figure [3j 



We can repeat the same procedure for the process in which the radiation vertex is inside 
the medium (see figure H]). For one scattering, 

cyg _ f d4 P A d 4 p! 4 d 4 qi 4 d 4 ki 4 ixi .( pi - p ) iaWq-qi) JVljk-h) iwfai+fci-pi) 

1 " J (2vr)4 d Xl (2vr) 4 V (2vr)4 d Xl (2vr)4 d 2/16 e 6 6 

X •"(g)^4A'A , 1 K+^ x l±) 2 } ■ ^'ij'" 2 ^4aiA(^1+)Xi±) 2 . i\40) 
1 or 4- -fii-fti rif 4. «£ o z 4- 



+ is ^i^ 11 + ^e p z + is 

xr;; ( A-)^/ aai V'^ , '(-fc,fci,o)4( 1 / 1+ ,y 1± )" / -' / " 1/ " 



k\ + ie 

(16) 

The gluon vertex 0, together with the metric from the propagator and the polarization 
vector can be simplified using k% + = k + , and so, the Dirac structure takes the form: 

«(?)^/(%^ (17) 

(b) (a) 

where 

(a) = 7_pi + ^4 1 p+7_ = pi + Aip_7„7 + 7_ = 2p 1+ A x p + ^^ , (18) 
(6) = 2 gi • - fi(g)^ 4[ ■ (19) 

~o 

In this last simplification, we are keeping only the dominant term of the Dirac equation 
since gi+ = g+. But we must not forget that the qj_ coming from the spinor u{q) is actually 
qix 7^ q_i_. This means that in the squared modulus, the transverse momentum that 
appears from this T-matrix corresponds to an inner momentum. The same is applied to the 
gluon transverse momentum coming from the gluon polarization vector, = — qu_. 

Using the properties listed in appendix 1X1 and 

dq- " < l-L_ e -i Q _(a; (i+1)+ -a; i+ )+iqx-(x (i+1)x -x i _ L )_ 1 



2n (2vr)2 g_ - (ql/2q + - te) (20) 



where 



V^ikxM, *s) = 9 ap (ki - fca) + g pd (k 2 - k 3 ) a + g da {k 3 - k 2 f 
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n l I \ Q+ j ip+ ( x (i+i) ~ x i)i 

Go(^ + ,x^;, (m)+ ,x (muk+ ) = 2m(x{m) _ Xj)+ exp|— {x{i+i) _ Xi) + 

/•ri(i(i + i) + )=Jt( j+ i)i rx{ i+ i)+ fdr ± \ 2 \ 



(21) 



is the Green's function of a free particle that propagates in the transverse plane from Xjj_ 
at (light-cone) time Xi + to X(j +1 )j_ at time we get for the T-matrix for a gluon 

emitted inside the medium 

' J 9 = \j d y+d*±e~ lCi± ' x± G A , A < i (y + ,y ± = ± ;L+,x. ± \q+)igT^ Ai 

x W Al A(x 0+ ,y+; ± )W aai (y+, L+- ± )u{q)f{k) 7 _M h {q + k), 

where the use of uppercase (lowercase) color indices in the Wilson lines indicate that they 
are to be taken in the fundamental (adjoint) as they correspond to the rescattering of a 
quark (gluon). 

The total T-matrix, 7 tot, is the sum of both contributions (equations (|15p and (|22p ). 
The spectrum is computed as the inelastic cross-section over the elastic cross-section (see 
the elastic process in figure [5]). Thus 



T tot | 2 



\M tot \A = X 7 = (|M||) + (|M||) + 2Re({M s il4}) , (23) 
where 

7 el = u(p) M h (p) = V2 P+ 1 M h (p)\ 2 . (24) 

As a consistency check we are able to recover the vacuum contribution in the limit of 
x — >• 1 from the quark amplitude, 

^_^_.)/i±CL=£>!} (25) 



dl 

x- 



dxd 2 \i 



± 



X— 5-1 



CfQ-s 1 _ a s 1 , - 

2vr 2 k ± 2 " 2vr2 k± 2 9 ^ q[X j 1 j 



with k^ = — qj_ and the vacuum splitting function 



P™ c q (z) = C F 



1 + (1 - S)2 " ^C F . (27) 



X— >1 



As for the other two terms in equation f|23|) (the medium contribution), the Dirac and 
color algebra are still to be simplified. They can be simplified using the polarization sum 
(with 77 = (0,1, ± )) 



£ e;(k, X)e v (k, A) = + k ^+^ (28) 



Figure 5: Elastic process. 
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and the relation between the Dirac spinors 

u (q, s) a u{q, s)p = f aj3 + m a/3 . (29) 

s 

Using these two relations, we will end up with the trace of 7- matrices that are easily 
computed. For the color algebra, one can reduce all the traces to the fundamental repre- 
sentation using [¥T] 

W ab ( X± ) = 2Tr [T a W F (x ± )T b W F Hx ± )\ (30) 
to simplify the expression. 

Putting all the kinematics in terms of the initial energy p + and the fraction of momentum 
carried away by the gluon, x, the medium amplitude can be written as: 



tmedl 2 ) = (\M g \ 2 ) + 2Re(Af fl Af} / 

= g 2 C F 1 + (1 ~ x)2 — Re(- \ f dy + dy + dx±dx ± dz ± e- ici± <^-*^ 

x p + y (1 - x)xp + J 

1 d 

x iv5y7 Tr ( G(y+,y± = °^y+> z ±\i+) w ^y+'y+-> ^) F 

Id 1 

Tr(G\y + ,y± = 0±; L + ,x ± \q + )G(y + ,z ± ; L + ;x_ L \q + ] 



Tiro— \ wt J -J- _l_ ? ~t~ 1 *-J- ) " \o-t- 1 "J- 1 *-T 1 -L "IT/ / 

TV ayj_ \ if 
+ 2 % • J dy+dx ± e^-- x± 
1 9 

x 



N8y ± 



Tr(G(y + ,y ± = Or,L + ,x ± \q + )WHy+X + ;0 ± ))\, 



(31) 



where the internal momenta were substituted by 



d 

qi±Go(y + ,y± = 0+ ±;x' 1+ ;x 1± \q + ) = i-^-Go(y + ,y ± = ± ; x 1+ ; x 1± \q + ). (32) 



Making the medium averages whose results can be taken from |42j . the medium ampli- 
tude for one gluon emission in the limit x — > 1 is, 

^P g <- q {x) Re{ — -L f dy + dy + dx ± e-^e-^fy L + +d ^^ 

P+ I (1 - x)xp + J 

ox± dy± 

+ 2^^ dy + d X± e- t(l ^-^-X(y + ,y ± =0 ± ;L + ,x±\q+)\ . 

(33) 



L med 

I I 

d d 



Although this expression was derived in the limiting case when x — > 1, we end up with 
an expression that depends explicitly on the fraction of energy carried by the emitted 
gluon. The result that we found for the medium amplitude is the exact vacuum splitting 
function, that appears as an overall pre-factor, corrected by two terms that depends on 
the medium parameters. Also, we have an explicit dependency with the initial energy that 
did not appeared in the limit of x — > 0. Moreover, if in this expression we take this limit, 
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x — > 0, exchange the momenta of the path integral, q+ by k + , exchange the transverse 
momenta by k^, and take the color representation in the adjoint one, we recover the 
previous results derived in the limit x — > (see ref. [2Z])- This indicates that it is possible 
to generalize this expression to include both limiting cases of the single gluon emission 
spectrum inside a medium. We will see a possible generalization in the next section. 



3 Multiple soft scattering approximation 

Being the general analytical solution of the path-integral from equation (j33[) unknown, an 
approximation scheme must be usecH. One of the possible choices is the dipole approx- 
imation. This is based in the observation that the Yukawa potential ao from the dipole 
cross section (eq. ([8])) has a leading quadratic dependency for small transverse distances 
|rj_| . This allows to write the dipole cross section to logarithmic accuracy as |13yi5| 

n(£>(r±)^«?(£)ri, (34) 

where q(£) is the transport coefficient which characterizes the transverse momentum 
squared , fi 2 , transferred from the medium to the projectile per mean free path, A. This 
is the main parameter of this approximation and encodes all the dynamical properties of 
the medium. Generally, this parameter is time-dependent since the medium is expanding 
but, for a static medium, it can be written as q = /U 2 /A. 

In this paper, we study only the case of a static medium, since the results can be easily 
generalized to an expanding medium^. In this case, the path-integral becomes that of a 
two-dimensional harmonic oscillator 1151. 



K(y+,y_L = 0±;y + ,x±\q+) = % OS c(y+,y± = 0±;y+,x±\q+) 

= — exp [iAi£i(x;[ +yi) - 2iA 1 x ± • yj , 

M = . lo Q ( +n V] , B X = cos My + - y + )\ , = ^/^, (36) 

2 sin [f%+ - y + )\ 2 ]j q + 

with imaginary frequency fL 

Note the subscript F in the medium transport parameter to distinguish it from the 
usual transport parameter that is obtained from the medium averages in the adjoint rep- 
resentation, 



L_Tr^( X± )^tA (y±) ^ =exp |_^i J dx + a(x ± -y ± Hx + )| 



N 2 

compared to the fundamental that is our case 



^Tr<V(x ± )W% ± )) =exp|-^ J dx + a(x ± - y± )n(x+)| . 



(37) 



(38) 



6 A comparison of the results for the soft gluon spectrum obtained using different analytical approxima- 
tions and exact numerical solutions for the path integral, can be found in [43]. 

7 The case of an expanding medium is just a generalization of equations (|3SJ and (f3"o]) (see refs. [271f4"4"] V 
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The two medium parameters are related by 



qF = c- A q ' 



(39) 



Substituting expressions ([35]) and (|36|) in (f33]) . one obtains the medium spectrum within 
the dipole approximation for the limiting case x — >• 1: 



dl 



with 



{Mmed) _ a s 

dq + dq ± 2 ~ 4(2vr) 2 ~ 4^ 
Re / dy+dy+dx.±e 



Pg^ q (x)(h + J 2 ) 



1 



x(l - 


- x)p+ 




9 r 




dy±~ 




1 



-2,4?Z) 



+ 



MfBi 2 



(D-iAxBi) 2 2(D-iA 1 B 1 ) 3 



q± 



exp 



-qi 



4(L> - MiBi 



representing the gluon term (< |M S | 2 )) and 



J 2 = — ^ . f dy + dx ± e- l ^^^-X(y + ,y ± = ± ;L+,x ± \E+) 



P+ Uo 



-Re< I dy + —^exp 



b: 



AA 2 B 2 



(40) 



(41) 



(42) 



the interference term. The variables Ai, B\ and f2 are the same as in equation (|36p . while 
A 2 , B 2 and D read 

A ^ 2 Bi n\a(f + -y +) ) ■ Bl=cosML + -«+)l , D = \(L + - (43) 



3.1 Interpolation between the hard (x — > 1) and soft (x — >■ 0) gluon 
emission 

To provide a heuristic interpolation between both limits, we should decide which variables 
we want for analyzing the spectrum since we have available q + , p + and x, but they are 
constrained by the relation q + = (1 — x)p + . The most direct single inclusive spectrum 
that we could compute from expression (|33p . should be the one related to the final quark, 
as previously derived. But since we want to find an interpolation function that also 
describes the single inclusive spectrum for the gluon, a good observable is the energy loss, 
independent of which parton in the splitting carries a larger fraction of the initial energy. 
Thus, if we substitute 

q + = (1 — x)p + — > x(l — x)p + = E + , (44) 
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when x — > 0, E + 
that 

cU 

q+ dq + dq A 



xp+ = k + and, when x — > 1, E + — > (1 — x)p+ = q. 

dl , dl 



This means 



d[x(l — x)p + ]d(\^ 



<l+ 



dk + dq± 2 

dl 
dq + dq± 2 



for x — > 0. 



for x —7- 1. 



(45) 



The transverse momenta are constrained by 

q± 2 < 2x 2 (l - xfp\ . 

As for the medium parameter, qp, we can substitute it by [14"ll21~] 

„ , \ ~ xC F + (1 - x)C A . 
g = xq F + (1 - xjgA = q A - 



(46) 



(47) 



With all these considerations, we are able to describe both limits within the considered 
approximations and recover the results that were derived in this manuscript and in [27j . 
So, applying the above formulas to equation (|40p we get 

dI med 



a s CF 

dE + dc^_ Ait 



with 



E 4 



h = ^2" Re / dy + dy + dx±e 



X{1 - x)P g ^ q { X ){h + I 2 ) 



(48) 



-"U-xj.p-2 Jy+ 



d d 

x — K(y + ,y_L = 0±;y+x±\Ej 



9x_l dyj 
—77- x Re 
-2.4fL> 



dy A 



(49) 



and 



(D-iAxBx) 2 2{D-iA 1 B x f 

dy±_ 



exp 



-qi 



4(L> - iAi-Bi) 



i* 8 



^ "7^2 eX P 
^2 



(50) 



4:A 2 B 2 

The variables Ai, A2, B\, B 2 , O and -D are the same as in equation (|43p but with the 
substitutions flUD and (|47|) : 



A 

A 2 = 



£+0 



2 sin - y + )\ 

E+n 



Bi = cos - y+)] 



f2 



1 - i 



1 



,, • ro , r TT > -Bi = cos [^(L+ - y+)] , D = -{L A 

2sin[S2(L+ - y + )\ 4 



>J_ 

E + 



(51) 



Note that while for either x — > or x — > 1, the meaning of qj_ in Eq. (|45p is clear 
as the transverse momentum of the emitted soft gluon (where q^ = — kj_) or quark 
respectively, the interpretation for intermediate x is far more involved. Actually, such 
kinematical situation requires the computation of a double inclusive cross section in which 
the Brownian motion of both outgoing partons is considered. We leave this computation 
for future studies. 
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4 Numerical Results 



The results for the double-differential medium-induced gluon radiation spectrum for a 
quark traversing a static medium in the limit of hard gluon emission are presented in 
figure The results are plotted as functions of the following dimensionless variables 
(sec [24,25] for the corresponding variables in the case of soft gluon emissions): 

ut = \<IaL\ , = ^72" , k 2 = -^j- , x . (52) 

2 ujJ qAL\ qAL+ 

The left plots correspond to the double differential spectrum. By integrating out qj_, we 
obtain the right-handed ones (the analytical expressions for this integration in the limit 



u c L + — > oo are written in appendix [B]) . Figure 6(a) presents the results for small medium 
parameters and figure [6(bJ| for larger ones. 

The two asymptotic behaviors that one expects are confirmed: when the energy or the 
transverse momentum of the emitted parton increase, the spectrum goes to zero. Another 
feature of the spectrum is the suppression in the low k? region. For values p+/u c ~ 1 
for every value of x, the spectrum becomes negative - the total radiation is smaller than 
in vacuum -, which does not happens when p^/ujf is smaller. For the case of the soft 
gluon emission, this suppression can be understood in terms of the formation time of the 
outgoing gluon: 

t9 form — ~T~ ■ ( 53 ) 



In the present case, we can generalize this concept to include both final particles and 
define the formation time of the outgoing parton as 

- 2E+ 

In the regime where tf orm 3> L+, the function becomes highly oscillatory and we enter 
in the regime of the Landau-Pomeranchuk-Migdal (LPM) effect [15,45]. As a result, a 
suppression of the spectrum due to the inefficiency of exchanging information with the 
medium is observed. In the upper left plot this behavior is less striking since we increase 
the medium length by two orders of magnitude. 



5 Conclusions and Outlook 



Summarizing, we have computed the medium-induced gluon radiation spectrum in the 
limit of a hard emission. Making use of the path-integral formalism to describe the prop- 
agation of high-energy particles inside a medium in terms of Wilson lines and Green's 
functions, we found a finite-energy correction to the double differential spectrum of the 
final quark. Since we are able to recover the soft gluon emission limit with the expres- 
sion derived in this manuscript, we provide a generalization that is able to describe both 
limiting cases. The result for this interpolation function (shown in equation (|48p for the 
multiple soft scattering approximation) contains, as in the case of soft gluon emissions, 
two contributions for the medium spectrum: the gluon term, which is suppressed by the 
initial momenta, and the interference term. The spectrum vanishes, as expected, in the 
kinematical limit. 
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K 2 




(a) 




Figure 6: Evolution of the medium- induced gluon radiation spectrum (|48|) of a quark in a 
static medium for different values of x and p + / w+ . The left plots correspond to the double 
differential spectrum (two different values of p + /uj^ in (a) and (b)) and the right ones are 
the same integrated over (two different values of w+L + in (a) and (b)). 
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Applying the multiple soft scattering approximation, we obtained explicit expressions 
for the case of a static medium. The integrals in the longitudinal variables were per- 
formed numerically and the main results were shown in figure [6) The conclusions are a 
clear suppression of the spectrum for small values of the parton transverse momentum, 
as predicted by the LPM effect. Also, the density of the medium constrains strongly the 
possible energy range of radiation that is emitted. 

Finite-energy corrections are a key ingredient for the quantitative description of the en- 
ergy loss processes in whatever used formalism. They are automatically included in Monte 
Carlo approaches, though until now in a heuristic form, see e.g. [28H30| for a discussion 
of its implementation and impact in a Monte Carlo framework. From this work we found 
a generalization that is still half-way of its true form: although we have an expression 
that is valid for all fraction of energies of the emitted parton, there is still one parton (the 
hardest one) that has its movement constrained in the transverse plane (this one is always 
described by a Wilson line). The outlook of this work is the implementation of these 
corrections in a Monte Carlo code in order to study its phenomenological consequences 
since now we are able to access the intermediate region of x, and, on the other hand, the 
full generalization of this result through allowing both final particles to have a Brownian 
motion in the transverse plane. 
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A Gamma matrices 

In the light-cone the gamma matrices read 




(55d) 



(55b) 



(55a) 



(55c) 



They obey to the following relations: 



( 7 +)t= 7 - ) 7_L'7_L = -2, 
( 7 -)t = 7 + 5 7 + 7 + = 0) 

(7_0 t = -7_L, 7~7~ = 0, 



7 7 7 = 27' 
7~7 + 7~ = 27" 
7_i_7"7_i_ = 27 



7_i_7 7_L = 27 + 
7 + 7_L7 + = , 
7~7_l7" = 0. 



(56) 
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B BDMPS limit 



From expression (|48|) . we are able to compute the BDMPS limit, taking the limit R = 
uj c L + — > oo. This is equivalent to performing the integration in the transverse momentum 
taking into account an opening angle for the parton emission, 0, < q± < xE+, where 
X = sinO and, then, take the limit x ~~ ^ °°- Doing the integration one finds that 



dl 



med 



dE 4 



4tt 



(l-X)xPg^ q {x){h+I 2 ), 



where 



h 



El 



- 1 + cxp 



HD-iAxBi) 



iA 1 B 1 ( X E + ) 2 
4(L> - iAiB^ 2 



SA 2 
B 2 E + 



exp 



AA 2 B 2 



1 



(57) 

(58) 
(59) 



and A\, A 2 , B\, B 2 , Q and D are defined in equation ([5T]) . 

Taking the limit x ~ * 00 ) the previous equations can be reduced to 



lim h = 2Re 



L+ , Qcos(r2yj 



o 



lim I 2 

R-¥oo 



-2Re 



dy^ 



sin(ily_| 

n 



(60) 



(61) 



cos(riy_)_) sin(f2y + ) 

Both integrals are divergent but, when summing the two contributions, the integral over 
the remaining longitudinal coordinate is finite and the the result is given by: 



dl 



med 



dE, 



2tt" 



x(l — x)P 9< _g(x)Reln [cos(fiL + )] . 



(62) 



This result agrees with the ones previously derived when doing the limit x — > 1 and x — > 
respectively. 
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